Connection between vibrational dynamics and topological order in simple glasses 
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We compare vibrational dynamics in two structurally distinct, simple monatomic model glasses 
simulated by molecular dynamics: the Lennard- Jones glass with an fcc-related structure and a glass 
with predominantly icosahedral short-range order. The former, characterized by a single local quasi- 
periodicity, supports only modes with acoustic behaviour. In the latter, the presence of optic modes 
and two incommensurate length scales is observed. This pattern of vibrational dynamics is shown 
to be closely related to that of a Frank-Kasper crystal having the same local topological order. 



PACS numbers: 63.50.+X, 61.43.Bn, 81.05.Kf 



I. INTRODUCTION 

The nature of vibrational dynamics in topologically 
disordered solids remains an unsolved problem of fun- 
damental interest. An important question is to what 
extent is the vibrational behaviour determined by the 
atomic structure in such materials. Liquids and glasses 
lack global symmetry, and, therefore, their vibrational 
eigenstates are no longer phonons. Nevertheless, long- 
wavelength phonon-like propagatingjacoustic plane waves 
can be observed in these structuresoH up to a certain fre- 
quency, known as the loffe-Regel crossover luir, beyond 
which they are heavily damped due to strong scatter- 
ing. However, the spectral functions (e.g. the dynamical 
structure factor S{Q,ui\\,&re still peak-shaped even be- 
yond u!iR (see e.g. RefJd'0 and references therein). Thus, 
the concept of dispersion, as peak-position frequency ver- 
sus the waygMajjtcyyiiagnitude, is still valid in disordered 
structuread'BErala'B'Ll. In glasses, it has been assumed that 
the position of the main maximum, Qq, of the static 
structure factor S{Q) plays the role of the first reciprocal- 
lattice point, wit|h|±he first pseudo-Brillouin zone extend- 
ing up to Qo/2m. 

However, the main maximum of S{Q), manifesting the 
dominant pseudo-periodicity in glasses, does not contain 
full information about the glass structure. Even in sim- 
ple monatomic liquids and glasses, there exist distinct 
types of short-range order (SRO) closely related to the 
respective counterpart periodic structureal3. Then, to 
what extent is the structural diversity of glasses reflected 
in their vibrational properties? In particular, can the 
distinctive features of the vibrational dynamics charac- 
teristic of a crystalline structure be discerned in the vi- 
brational properties of a glass with similar local order? 

A type of SRO that is most interesting in this respect 
is that based on icosahedral coordination of the first shell 
of neighbors. It is commonly associatedi-with the glass- 
forming ability of simple metallic systemstil. The frustra- 
tion inherent in packing icosahedra in Euclidean (flat) 3D 
space in the liquid state produces 6-fold defects which. 



on cooling, form disclination lines. The latter may align 
in space, forming globally ordered patterns, e.g. Frank- 
Kasper periodic phases or morphologically related qua- 
sicrystals. Under sufficiently rapid cooling, the disclina- 
tion lines get entangled, forming a glassM. 

The large numbers of atoms in the primitive cells of the 
Frank-Kasper phases give rise to many optic vibrational 
modes. Can these modes survive the lack of long-range 
order in the related glass? In addressing this issue, it 
is important to discriminate between the optic modes 
arising as a result of a particular type of topological SRO, 
of interest here, and those induced by the chemical SRO 
in a system with more than one type of at(MHp-|The latter 
have been observed in simple glassy alloyalSEj. 

Here, we employ molecular dynamics to investigate the 
vibrational dynamics in a monatomic model of a glass 
with icosahedral SRO (IC glass)Eil. The thermodynam- 



ically stable solid phase 
dodecagonal quasicrys 
structure is the a phasi 



freezing for this system is a 
hereas the energy- favoired 
13, a Frank-Kasper crystaO. 
The rest of the paper is arranged as follows. In Sec. II, 
we describe simulation procedure. The results and their 
discussion are presented in Sec. [II and Sec. respec- 
tively. The conclusions are given in Sec. 0. 



II. SIMULATION 

The IC glass was obtained in the course of a NVE 
molecular dynamics (MD) simulation. A system of 16000 
particles was cooled from the liquid state above the melt- 
ing point Tm — 0.5 at a constant density p = 0.8771 (all 
quantities reportedphere are expressed in the Lennard- 
Jones reduced unitsEJ defining the pair potentialifl) . This 
is the rninimiim energy density at zero pressure in the 
a phasell3. We cooled the system in a stepwise man- 
ner. At each cooling step, the temperature of the sys- 
tem was reduced by scaling the velocities of the con- 
stituent particles. Then the system was equilibrated at 
constant temperature until the decrease of the poten- 
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TABLE I: Temperature dependence of the statistics of —72° 
disclination lines in the IC glass. 



FIG. 1: (a) Frank-Kasper ZlA polyhedron, (b) A fragment 
of the —72° disclination line (two interpenetrating ZlA ele- 
ments) as observed in the simulated IC glass. 



tial energy was observed to cease, indicating that a state 
of metastable equilibrium had been reached. Successful 
equilibration of this system in the metastable Uquid do- 
main below the melting point was possible due to the 
complexity of its crystallization pattern, significantly in- 
hibiting the process of crystalline nucleation. The equi- 
libration time rapidly increased with decreasing temper- 
ature. At T ~ 0.3, this time exceeded the time of the 
simulation run Below T = 0.3, we found that the equili- 
bration had not been completed after the relaxation run 
of t = 10'' At, At = 0.01. At this point, the system, 
therefore, was assumed to be in the glassy state. Then 
it was cooled to the temperature T — 0.01 at which the 
atoms vibrate around their equilibrium positions and do 
not diffuse away, and their vibrational dynamics was in- 
vestigated. 

Mindful of the possible growth of crystalline grains in 
the process of glass formation, we analyzed the evolu- 
tion of the six- fold structural elements, the Frank-Kasper 
Zli polyhedra (shown in Fig. 0(a)). In the process 
of Frank-Kasper crystallization, interpenetrating ZlA el- 
ements form —72° disclination lines (Fig. |l|(b)) which 
eventually thread the entire samplalil. Therefore, the ex- 
tent of crystallization can be assessed by monitoring the 
length of these disclination lines. Table | demonstrates 
that the process of glass formation in the IC system is not 
accompanied by any significant growth of the crystalline 
phase. 

As a reference, we simulated a monatomic Lennard- 
Jones (LJ) glass with the SRO related to the face- 
centered cubic (fee) structure. This glass was prepared 
at the density p = 1.0888, corresponding to the zero- 
pressure energy minimum of the fee crystal with the LJ 
potential truncated at Tc — 4, by rapidly quenching a 
system of 16000 particles from a well-equilibrated liquid 
state at T = 3 to T = 0.01. The structural distinction 
between the two glasses is demonstrated in Fig. ||. 

In order to gain an insight into the vibrational prop- 
erties of the IC glass, we also studied, the dynamics of a 
one-component model of the a phaseEO. The latter model 
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FIG. 2: Static structure factors S{Q) for the IC glass (solid 
line) and the LJ glass (dashed line). The inset shows the 
small-Q part of S{Q) (solid line, left-hand scale) and the par- 
tial structure factor calculated for the centers of icosahedra 
(dash-dotted line, right-hand scale), both for the IC glass. 
{QoP^^^'^ — 7.17 and 7.08 for IC and LJ glass, respectively.) 



utilized the same pair potential as the IC-glass model. 
Both models were studied at the same density. 

We associated longitudinal vibrational modes 
with non-zero frequency maxima in the dynami- 
cal structure factor, S{Q,lu), the spectrum of 
density-fluctuation time-autocorrelation functl 



3f tiie 
tionkil: 



S{Q,uj) = N-^J^dt e''^*(p(Q,t)p(-Q,0)), where 

^(Q)*) = EsLi'3^P[~*Q ■ ^k{t)] and rfc(t) is the 
positional vector of the particle k. For transverse 
vibrational modes, we calculatedi— ithe spectra of 
the transverse current fiuctuationa22l: Ct{Q,i-u) = 
Q'^N-^j^dt e'"*(it(Q,i)jt(-Q,0)), where jt(Q,t) is 
the transverse component of the Fourier-transform of the 
local current: jt{Q„t) = X]f=i(et •Vfc(t)) exp[-zQ-rfc(t)], 
Gf _L Q and Vfc(t) is the velocity of particle k. Finite-time 
truncation effects in the spectra of the time-correlation 
functions were reduced by convolution with a Gaussian 
window function. 
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FIG. 3: Dynamical structure factors in the simulated systems at three values of wavenumber as marked, (a)-(c): IC glass, solid 
line; dashed line, a phase. In (b), the dotted curves show two Gaussians, the sum of which is fitted to the solid curve; this fit, 
if plotted, is indistinguishable from the solid curve, (d)-(f): LJ glass. 




FIG. 4: (a) Dispersion of vibrational excitations for the LJ glass. Open squares: longitudinal branch. Open circles: transverse 
branch. Dash-dotted line: the first Brillouin-zone boundary for the fee structure in the (111) symmetry direction, (b) Vibra- 
tional density of states for the LJ glass (solid curve) and for the fee structure (dashed curve), (c) Phonon-dispersion curves for 
the fee structure in the (111) symmetry direction. 



III. RESULTS 

Fig. ^ shows S{Q, ljj) for the two glassy models, as well 
as the a phase, for several representative values of Q. For 
sufficiently small Q, both glassy systems clearly exhibit 
acoustic modes which manifest themselves as pronounced 
Brillouin peaks in S{Q,uj), Fig. |^(a). A marked differ- 



ence between the two sets of curves arises for Q > 0.2 Qq. 
The width of the Brillouin peak in the LJ glass becomes 
comparable with that of the whole vibrational spectrum 
(see Fig. ||(e)) and for larger Q the peak disappears 
(Fig. ^(f)). For the IC glass, a well-defined Brillouin 
peak persists well beyond Qo (Fig. ^(c)). Another im- 
portant feature of the IC glass is that, for Q > 0.2 Qq, 
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FIG. 5: (a) Dispersion of vibrational excitations in the IC glass. Open squares: longitudinal branch. Open circles: transverse 
branch. Solid circles and triangles: maxima positions of the two Gaussians approximating S{Q,ui). Dash-dotted line: the first 
Brillouin zone boundary for the cr-phase structure in the (001) symmetry direction, (b) Vibrational densities of states of the 
IC glass (solid line), and the cr-phase structure (dashed line), (c) Phonon-dispersion curves for the a-phase structure in the 
(001) symmetry direction. 



the Brillouin maximum develops a well-defined shoulder 
on the high-frequency side (Fig. ^b)). We analyzed this 
structure of the Brillouin peak by fitting it with the sum 
of two Gaussians. 

The assumed interpretation of the shoulder observed in 
the Brillouin peak in S{Q,uj) of the IC glass as indicating 
the presence of a separate maximum can be corroborated 
by a comparison with the structure of S{Q,uA calcu- 
lated for the a phase which is also shown in Fig. ^(a)-(c). 
S{Q, Lu) for the a phase was computed at T = 0.8 in the 
same way as for the glassy systems: the averaging was 
performed over all wavevectors Q confined to a narrow 
spherical shell in reciprocal space, such that |Q| ~ Q. For 
Q/Qo — 0.35, S{Q, iu) exhibits a pronounced peak at the 
frequency close to the position of the shoulder observed 
in S(Q, oj) of the glassy phase. The two sets of curves in 
Fig |3|(a)-(c) clearly show certain similarities: the number 
of peaks is the same and their corresponding positions 
are close. This leads us to the conclusion that the shape 
of the dynamical structure factor in the frequency do- 
main for the IC glass is generically related to that for 
its crystalline counterpart, the a phase. Namely, the 
peaks in S{Q,uj) for the IC glass originate from the cor- 
responding peaks for the a phase, which are broadened 
and shifted by disorder-induced level-repelling effectsE3. 
The level-repelling and hybridization effects describe a 
reconstruction of the crystalline vibrational spectrum to 
give the disordered one. These should give rise to a shift 
of the low-frequency crystalline peaks downwards in fre- 
queuCiy but upwards for the high-frequency peaks (see 
Ref.B for more detail). Such tendencies are indeed seen 
in Fig. 11(a)- (c) for the relative positions of the peaks for 
the a phase and the IC glass. Therefore, our results sup- 
port the conjecture that disorder-induced level-repelling 
effects in vibrational states in crystalline counterparts of 



topologically disordered structures, such as the IC glass, 
are responsible for the main features of their vibrational 
spectrum. 

The dispersion relations as derived from the positions 
of the respective Brillouin maxima for the LJ and IC 
models are shown in Figs. ^ and ^, respectively. The 
patterns corresponding to the longitudinal modes in the 
two glassy systems are strikingly different. The longi- 
tudinal dispersion curve for the LJ glass demonstrates 
a pattern typical for acoustic modes, with a maximum 
centered at around Qa/2, Fig. ^a). The dispersion pat- 
tern of the single Brillouin peak observed in S{Q,ll!) for 
the IC system for Q < 0.2 Qo can also be identified with 
acoustic behavior. For larger values of Q, we identify two 
separate dispersion branches associated with the posi- 
tions of the maxima of the two Gaussians approximating 
S{Q, uj), as shown in Fig. ||(a). The main subpeak of the 
Brillouin maxima exhibits a very weak Q-dependence, re- 
sembling the behavior of the optic modes in crystalline 
solids (Fig. 1^). The shoulder observed in the Brillouin 
peaks of the IC system gives rise to a higher-frequency 
dispersion branch which exhibits a maximum. However, 
in contrast to the acoustic-type dispersion observed in 
the LJ glass, the maximum is located at around 0.4 Qo- 



IV. DISCUSSION 

The pronounced differences in the behavior of the lon- 
gitudinal modes in these two model glasses can be ratio- 
nalized by considering the respective modes in the crys- 
talline counterparts - the a phase in the case of the IC 
glass, and the fee crystal in the case of the LJ glass. 

The LJ glass is characterized by a single dominant 
structural repeat distance, manifested by the first peak 
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of S{Q). It corresponds to the interplanar spacing along 
the (111) symmetry direction in the fee structure (i.e., 
the apex-basal-plane distance in a close-packed tetrahe- 
dron). Qq being the first reciprocal-lattice point, the 
first Brillouin-zonc boundary is at Qo/2, which corre- 
sponds to the maximum in the dispersion curve of the 
longitudinal acoustic (LA) modes (Fig. |4|(a)). As a re- 
sult of this structural simplicity, the vibrational disper- 
sion pattern in the LJ glass is also simple, with only LA 
and transverse branches, as in the counterpart crystalline 
(fee) structure (Fig. |4|(c)) that has just one atom as the 
basis in the primitive unit cell. 

The pattern of longitudinal vibrational modes in the IC 
glass is considerably more complex (Fig. |^(a)), displaying 
a low-frequency weakly dispersive branch with the first 
maximum at Q/Qo — 0.2. Its dispersion behavior clearly 
resembles that of the optic modes which dominate this 
range of frequencies in the vibrational spectra of the a 
phase (Fig. |5|(c)). In the same domain of Q, the IC glass 
also exhibits a higher-frequency branch associated with 
the high-frequency shoulder of the Brillouin peak. As has 
been mentioned earlier, this shoulder mirrors a peak in 
S{Q, uj) of the high-temperature a phase (Fig. ||(b)). The 
latter manifests optic modes and therefore we conclude 
that optic-type modes dominate the vibrational dynam- 
ics in the IC glass for Q/Qo > 0.2. The position of the 
maximum in the high-frequency branch, Q/Qo — 0.4, 
is twice that in the low-frequency branch, Q/Qo ~ 0.2. 
However, both branches are suppressed at Q = Qo by 
the influence of the main peak of S{Q). 

The appearance of modes in the IC glass, mirroring 
the optic modes in the corresponding cr-phase crystal 
(Fig. 11(c)), is a consequence of the type of topological 
SRO that exists in both structures. The first Brillouin- 
zone boundary in the a phase in the (001) symmetry 
direction is at Q ~ 0.2 Qo (Fig. ||(c)), corresponding 
very closely to the first pseudo-Brillouin-zone boundary 
value found in the IC glass (Fig. |5|(a)). The behavior 
observed in the IC glass, that the low- frequency branch 
first peaks at Q ~ 0.2 Qo whereas the high-frequency 
branch peaks at Q ~ OAQo, is typical for nonsymmor- 
phic crystalline structures, of which the a phase is an 
example (space gro«p P42/mnm), for Q in the direc- 
tion of a screw axisE3. Some modes sense any particular 
(pseudo-) periodic atomic-density fluctuations, whereas 
some modes sense alternate such atomic-layer spacings 
and therefore the apparent first Brillouin-zone bound- 
aries differ by a factor of two. On the other hand, the 
dominant length scale in the IC glass that is manifested 
by the main peak of S{Q) corresponds to the layer pe- 
riodicity in the (111) direction of the a phase (the same 
apex-basal-plane distance in tetrahedra as in fee). 

Optic modes in the IC glass give no contribution to 
S{Q,uj) for Q/Qo < 0.2. This is not surprising be- 
cause the dynamical structure factor describes the re- 
sponse of the system to plane wave excitations. Even in 
crystals, the optic phonons are not simple plane waves 
(the atoms inside the unit cell are not displaced accord- 



ing to a plane wave), and in the long- wavelength limit 
{Q 0) i^gly acoustic modes respond to plane-wave 
excitationsoH. The signature of optic-like excitations ap- 
pears in the dynamic structure factor of the IC glass only 
when the wavelength becomes comparable with the size 
of structural units responsible for this effect. 

The dispersion curves for the IC glass thus indicate 
the existence of an additional pseudo-periodic interlayer 
spacing incommensurate with that manifested by the 
main peak of S{Q). The latter corresponds to the same 
pseudo-period in the atomic-density fluctuations as for 
the LJ glass. The second length scale is related to the 
interlayer spacing in aggregates of interpenetrating icosa- 
hedral units, as along the (OOl)-axis in the a phase. This 
assertion is confirmed by the fact that both the total 
S{Q) and the partial structure factor calculated only for 
the centers of icosahedra exhibit a prepeak at Q ~ 0.4 Qo 
(see the inset in Fig. ||). The presence of an additional 
length scale in the IC glass appears to agree with recent 
observations of med.itun- range order in paracrystalline 
amorphous modelsEJuJ. This length scale is associated 
with a particular type of structural unit which is mor- 
phologically close to those found in the counterpart crys- 
talline structure. In a reciprocal-space picture, the effect 
of these structural u]iiita--pan be analogous to that of the 
"quasi-Bragg planes"E30 used for explanation of the first 
sharp diffraction peak in S{Q) of amorphous solids. 

The difference in the behaviour of collective vibrational 
modes in the two glasses is also reflected in their densities 
of states, Figs. ||(b), ||(b). On the other hand, the curves 
for the vibrational density of states in both glasses appear 
to be consistent with the respective curves calculated for 
the crystalline counterparts. 

The dispersion patterns of the transverse modes are 
qualitatively similar (Figs, ^(a), ||(a)). This implies that 
either these modes are insensitive to the structure or, 
more likely, any structure effects cannot be shown using 
Ct{Q,oj). This will be analyzed elsewhere. 



V. CONCLUSION 

In conclusion, we have demonstrated that structural 
distinctions in simple monatomic glasses give rise to an 
appreciable difference in the behaviour of vibrational ex- 
citations. We have shown that the main features in the 
dynamical structure factor of the IC glass are generi- 
cally related to those of its crystalline counterpart, the a 
phase. The other distinctive features of vibrational dy- 
namics displayed by the IC glass, namely, the presence of 
optic-like modes and an additional pseudo-Brillouin zone 
boundary incommensurate with Qo, can be accounted for 
by the spcciflc structure of aggregates of interpenetrating 
icosahedra characteristic of the a phase, which was pre- 
viously showrO to be a good crystalline counterpart of 
the IC glass. This observation implies that the structural 
similarity between a glass and its crystalline counterpart 
extends well beyond the shell of nearest neighbors. The 
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existence of an additional incommensurate length scale 
and optic modes may be a generic feature of glasses hav- 
ing complex counterpart crystalline structures. 
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